Quantum Hall ferromagnets, cooperative transport anisotropy, 
and the random field Ising model 



(N 
O 

o 

wo: 
< 



i 

c 

o 
o 



> 

(N 

o 

oo 
O 
(N 
O 
-i— > 

a 



i 

C 

O 

o 



J.T. Chalker 1 - 2 , D. G. Polyakov 3 |] F. Evers 3 , A. D. Mirlin 1 ' 3 ^ and P. Woffle 1 ' 3 
1 Institut fur Theorie der Kondensierten Materie, Universitdt Karlsruhe, 76128 Karlsruhe, Germany 
2 Theoretical Physics, University of Oxford, 1 Keble Road, Oxford 0X1 3NP, United Kingdom 
3 Institut fur Nanotechnologie, Forschungszentrum Karlsruhe, Karlsruhe 76021, Germany 

(Dated: February 1, 2008) 

We discuss the behaviour of a quantum Hall system when two Landau levels with opposite spin 
and combined filling factor near unity are brought into energetic coincidence using an in-plane 
component of magnetic field. We focus on the interpretation of recent experiments under these 
conditions [Zeitler et al, Phys. Rev. Lett. 86, 866 (2001); Pan et al, Phys. Rev. B 64, 121305 
(2001)], in which a large resistance anisotropy develops at low temperatures. Modelling the systems 
involved as Ising quantum Hall ferromagnets, we suggest that this transport anisotropy reflects 
domain formation induced by a random field arising from isotropic sample surface roughness. 

PACS numbers: 73.43.Cd, 73.21.-b, 75.10.-b, 75.10.Nr 



Two very striking experimental observations of large 
electronic transport anisotropy for quantum Hall syste 
in tilted magnetic fields have been reported recently. 
In both cases, anisotropy appears at integer values of the 
filling factor v with an in-plane magnetic field compo- 
nent tuned to bring two Landau levels of opposite spin 
into energetic coincidence. While the in-plane magnetic 
field component itself defines an axis within the sample, 
the fact that large anisotropy appears in resistivity only 
below a characteristic temperature of about 1 kelvin sug- 
gests it has a cooperative origin. Our aim in this paper 
is to develop a theoretical treatment of such systems and 
to discuss the source of the observed anisotropy. 

In view of the phcnomcnological similarities, it is nat- 
ural to make comparisons between these Landau level 
coincidence experiments and the earlier discovery of re- 
sistance anisotropy in quantum Hall systems near half- 
filling of high Landau levels,cl attributed to the formation 
of a uniaxial chacffe density wave with a period set by the 
cyclotron radius.u Some distinctions are, however, clear. 
Most importantly, the nature of the electron states near 
the chemical potential and their average occupation is 
quite different in each case: two separate orbital Lan- 
dau levels with opposite spin and a combined filling fac- 
tor close to unity are involved in the coincidence exper- 
iments, as against a single, spin-polarised and roughly 
half-filled Landau level in the other case. In this context 
it is desirable to examine a range of possible explanations 
for low-temperature anisotropy. 

The study of cooperative effects in coincident Landau 
levels has a long history. Consider a system with fixed 
magnetic field strength B± perpendicular to the two- 
dimensional electron or hole g (IS j CIS 3b function of total 
field strength -Btot- In a single-particle description there 
are pairs of Landau levels having opposite spin orien- 
tations and orbital quantum numbers differing by AN, 
which are separated in energy by ANhiu c — <?*/iB-Btot, 
where uj c oc B± is the cyclotron frequency and g* ^B to t 
is the Zeeman contribution. This energy gap falls to zero 
at coincidence. Inclusion of exchange interactions leads, 



at a combined filling factor of unity for the crossing lev- 
els, to a first-order transition between two ground states 
in which one or other level is completely filled. Within 
a Hartree-Fock treatment, the excitation gap remains 
non-zero through this transition^ Early observations of 
a quantised plateau in Hall resistivity p xy and deep mini- 
mum in diagonal resistivity p xx , both persisting through 
the transition,El provide support for such a picture, while 
measurements of the temperature dependence oL.p xx al- 
low determination of an activation energy gapEl which 
has the variation with B\\ expected theoretically. This 
Landau level coincidence transition is one example of a 
broad class of cooperative phenomena in quantum Hall 
systems, involving ferromagnetism of either spin or pseu- 
dospiii^ariables, which have been a focus for much recent 
work.uEl From such a perspective, representing the two 
Landau levels involved using two states of a pseudospin, 
interactions between pseudospins are ferromagnetic with 
Ising anisotropy,tJ while Btot measured from its value at 
coincidence acts on pseudospins as a Zeeman field. 

With this background in mind, we return to a dis- 
cussion of transport anisotropy. Following the sugges- 
tions of Refs.[l] and 0, it is clear that the presence 
of a spin or charge density wave could potentially ex- 
plain this observation. From the original Hartree-Fock 
calculations which were for a one-band model and con- 
sidered only trial states with homogeneous charge den- 
sity, it was found that instability to a spin density wave 
is preempted by the first-order coincidence transition. 
An escape from this conclusion might be provided by 
the fact that the sample involved is in one experiment a 
Si/GeSi heterostructure,El and in the other a wide quan- 
tum welhcl the former has valley degeneracy and the lat- 
ter has two occupied subbands. Alternatively, it might 
be that a better trial state in Hartree-Fock theory, or 
calculations which go beyond this approximation, would 
yield a stripe phase as the true ground state near coin- 
cidence. Howeiier, recent Hartree-Fock calculations for 
bilayer systemsliil give only ferromagnetic pseudospin or- 
der with parameters relevant in the present context, as do 
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calculations for one-band models, in which Hartree-Fock 
solutions with bath spin and charge density modulations 
are consideredjla or the Hamiltonian with realistic in- 
teraction potentials as diagonalised exactly for a small 
number of electrons.EJ 

The observations of Refs.|l| and || therefore present 
a puzzle, which we argue in the following can be un- 
derstood in terms of domain formation, with a char- 
acteristic size much larger than the relevant scale for 
stripe phases, the cyclotron radius. Our account in- 
volves three distinct ingredients. First, we suggest that 
domains are induced by a random Zeeman field acting 
on the pseudospins, which arises from the interplay be- 
tween isotropic sample surface roughness and the in- 
plane component of magnetic field. Second, we show 
that a random field generated by this mechanism is in- 
trinsically endowed with anisotropic correlations, and 
that the correlation anisotropy is large enough to ex- 
plain the observed anisotropy in resistivity. Third, we 
argue that transport in a multi-domain sample occurs 
along domains walls, via the processes discussed recently 
in Refs. p^[l5|Jl6| . The onset temperature for transport 
anisotropy arising by this mechanism is the Curie temper- 
ature of the Ising oiAantum Hall ferromagnet, and we note 
that the reportedlffl onset temperature of about 1 Kis 
similar to the value for tiws Curie temperature expectedt£l 
and observed elsewhereO 

As a starting point for a more detailed discussion, con- 
sider an energy functional for the system. Introducing co- 
herent state creation operators c j (r) and c j (r) for the two 
Landau levels involved, correlations are characterised by 
the expectation value of pseudospin, S(r) = (c^er^c^), 
where er is the vector of Pauli matrices. The order pa- 
rameter has magnitude |S(r)| = S, where S = 1 at a 
combined filling factor of unity for the two Landau levels 
and is smaller otherwise. For variations of S(r) which are 
smooth on the scale of the cyclotron radius, one expects 
the energy functional to have the form 

E = J (-DS 2 + J\VS\ 2 + SJ\d n S(r)\ 2 ~hS z )d 2 r. 

(1) 

Here D > represents Ising anisotropy, J is the spin stiff- 
ness, the derivative d„ = n- V acts in the direction of the 
in-plane magnetic field, denoted by n, and 5 J represents 
spatial anisotropy in the spin stiffness (for simplicity, we 
omit anisotropy in spin-space from the stiffness). The 
effective Zeeman field acting on pseudospins is h. In ex- 
periment, the strength of this field varies through zero 
as the tilt angle of the sample in the applied magnetic 
field is varied through the Landau level coincidence point; 
its strength depends also on field magnitude i?tot and on 
carrier density n. 

For a homogeneous system, the ground state of Eq. (|l|) 
is uniform with S z = sgn(h)S and Sj_ = 0. Domains 
may arise either in metastable states or because they are 
induced by quenched disorder. While metastability and 
hysteresis are .observed in some examples of quantum Hall 
ferromagnetsJLH this is not reported to be an important 



aspect of observations in Refs. [l| and g. We therefore turn 
to domains induced by disorder. Potentially the most im- 
portant source of disorder in Eq. (|l]) is randomness in h, 
and behaviour of the random field Ising model has been 
studied very extensively.Ea It is useful to distinguish the 
weak and strong disorder regimes: taking h to fluctuate 
about mean value zero with amplitude A and correla- 
tion length I, and supposing I is greater than the domain 
wall width w — W J / 'D, the boundary between the two 
regimes lies at I A ~ \J JD. At weak disorder, domain size 
£ is much larger than I and domain morphology depends, 
amongst other things, on the difference in energy per unit 
length of domain walls running parallel or perpendicular 
to the in-plane magnetic field component: this energy dif- 
ference is of order V SJD. At strong disorder, which we 
shall argue is the limit relevant here, the domain pattern 
is simply that of sgn(/i). 

To apply these ideas, it is necessary to identify a mi- 
croscopic origin for such a random field. One possibility 
is that variations in carrier density n, arising either from 
impurity scattering or from large-scale inhomogeneities, 
produce changes in the value of h. Randomness of this 
kind is expected to be spatially isotropic, but may give 
rise to transport anisotropy via dependence of the domain 
wall energy on spatial orientation. A second possibility 
is that sample surface roughness changes the local value 
of 9, and hence h. To compare the likely importance of 
these two, we appeal to experiment, noting (for example, 
from Fig. 2 of Ref . [l]) that while the coincidence transition 
has a rather small width (0.5°) in 9, it has a much larger 
width (20%) in i?tot, which is indicative of its width in 
n. The existence of sample surface roughness of the re- 
quired amplitude, with I ~ 1/im, is noted in Ref.|l| and 
well-established in a variety of other contexts. Ej More- 
over, it can account for transport anisotropy, as we now 
show. 

Let z(r) denote height of the sample surface above an 
average reference plane, and let 9 C be the critical angle 
at which the Landau level coincidence transition occurs. 
Then for small-angle roughness 

h(r) = a(9 - 6 C ) + ad n z(r) , (2) 

where a is a proportionality constant. Crucially, by this 
mechanism surface roughness with a correlation function 

C(r) = (z(r')z(r + r')) - (z(r')) (z(r + r')) (3) 

which is isotropic generates a random field with spatially 
anisotropic correlations, since 

(h(r')h(r + r')) - (h(r')) (h(r + r')) = -d 2 n C(r) . (4) 

To establish the characteristic degree of this 
anisotropy, we have carried out numerical simulations. 
Taking z(r) to be a superposition of overlapping Gaus- 
sian functions of position, with centers placed at ran- 
dom points in the plane and amplitudes distributed uni- 
formly about zero, and setting 9 — 9 C , the resulting ran- 
dom field h(r) in a typical realisation is illustrated in 



FIG. 1: Greyscale plots of: (upper panel) spatially 
isotropic surface roughness z(r), and (lower panel) the spa- 
tially anisotropic Ising model random field d n z(r) which re- 
sults from this surface roughness. 

Fig. |l]. Anticipating our discussion of transport on do- 
main walls, we quantify the anisotropy evident in this fig- 
ure by following the classical dynamics of a particle that 
moves along, contours of h(r), using methods described 
previouslyEI Averaging over randomly-placed starting 
points, we expect diffusive motion in the sense that mean 
square displacements grow linearly in time. Taking n 
parallel to the y-axis, the quantities D xx (t) = (x 2 (t))/t 
and D yy (t) = (y 2 (t))/t should then approach the eigen- 
values, D xx and D yy , of the diffusion tensor, for times 
t which are large compared to the correlation time, t . 
Evidence that D xx (t) and D yy (t) indeed tend to a fi- 
nite limit, with D xx ~ 8D yy , is presented in Fig.|^. The 
orientation of this anisotropy, with the larger diffusion 
constant in the direction perpendicular to the in-plane 
magnetic field component, is as observed in Refs.[l],[| and 
its magnitude is about the same as that determined at 
low temperature using a Hall bar sample .EJ The precise 
value of D xx /D yy will be dependent on disorder distri- 
bution, but we see no reason to expect large variations. 
Our calculations also provide an opportunity to test the 
universality class of our anisotropic percolation prob- 
lem, since diffusive growth in mean square displacement 
arises at long times from a balance between bounded 
motion on closed trajectories and super-diffusive motion 
on trajectories which remain open up to the observation 
time.Ea Averaging only over open trajectories, one ex- 
pects (x 2 (t)) oc (y 2 (t)) oc £ 8 / 7 , if the anisotropic problem 
is in the same universality class as the standard classical 
percolation problem: data shown in the inset to Fig.^ 
support this conclusion. 

The foregoing discussion is based on the idea that 
transport occurs along boundaries between domains. In 
order to substantiate this, we next examine transport 
properties of domain walls between oppositely magne- 



FIG. 2: Simulation data used to determine diffusion coeffi- 
cient anisotropy. Mean square displacements per unit time 
in directions perpendicular (D xx (t), full line) and parallel 
(D yy (t), dashed line) to the in-plane field direction n, aver- 
aged over all trajectories. Inset: averages over open trajecto- 
ries only, (x 2 (t)) /t 8 ^ 7 and (y 2 (t)) /t 8 ^ 7 , demonstrating scaling 
with the classical percolation exponent value. 

tised phases of the Ising quantum Hall ferromagnet. Re- 
calling that the domain wall forms the boundary between 
a region on one side with filling factors for the coincident 
Landau levels of ff ~ 1 and ~ 0, and a region on 
the other with interchanged filling factors, v\ ~ and 
Ul ~ 1, the simplest structure one might imagine is that 
shown in Fig. |](a) ■ In this picture, the wall supports two 
counter-propagating modes with opposite spin polarisa- 
tion, which arise as edge states of the occupied Landau 
levels in the domains on either side. Such an Ising domain 
wall, in which Sj_(r) = everywhere and S(r) = at 
the wall centre, may be stabilised by short-papge scatter- 
ing, which allows solutions with |S(r)| < 1,E3 in contrast 
to Eq. ([j]) . For a sample without shptt-range scattering, 
however, Hartree-Fock theory yieldsta the Bloch domain 
wall structure shown in Fig.||(b). Here, Sj_(r) ^ within 
the wall. In consequence, within Hartree-Fock theory 
there is mixing and an avoided crossing of edge states 
arising from occupied Landau levels on either side of the 
wall. 

At this level of approximation, for a Bloch wall the 
chemical potential lies within a quasiparticle gap. To 
account for transport under these conditions, it is nec- 
essary to consider collective excitations. The combined 
consequences of continuous symmetry for the Hartree- 
Fock solution under rotations of (S_i_(r)) about the Ising 
axis, and the connection between spin or pseudospin and 
charge that is standard for quantum Hall ferromagnetsEl 
have been examined in a related context in Ref.|l5|. In- 
troducing pseudospin rotation angle ip, as a function of 
position coordinate x along the wall and imaginary time 
t, the action 



S = 



dip\ 2 1 / dip^ 



dx 



dx dr 



(5) 



is obtained for domain wall excitationsiij where in our 
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(a) 




FIG. 3: Schematic summary of domain wall structure, show- 
ing pseudospin and excitation energies E, as a function of 
position y across the wall, within Hartree-Fock theory: (a) 
for an Ising wall stabilised by short-range scattering, and (b) 
for a Bloch wall. 

context p ~ Jw ~ e 2 /e and u ~ e /eh. A charge density 
(2tt)^ 1 dip / dx is associated with these modes. This is the 



action for a spinless Luttinger liquid. A vital property 
for our argument is that left and right-moving excita- 
tions propagate independently, provided rotation sym- 
metry about the pseudospin easy axis is exact. In short, 
Fig. ^|(a) remains a useful picture even without short- 
range disorder to stabilise an Ising wall, provided only 
that there is no spin-orbit scattering. In this picture, 
transport in a multidomain sample occurs at domain 
boundaries, via two independent, counter-propagating 
sets of modes. Neglecting quantum interference effects, 
we arrive at the problem for which numerical results are 
given above. 

In conclusion, we have argued that the observations of 
anisotropic transport reported in Refs.|l| and ^|can plausi- 
bly be attributed to formation of anisotropically shaped 
domains, induced as a result of sample surface rough- 
ness. Our numerical work demonstrates that this mech- 
anism generates, an anisotropy comparable to that found 
experimentally!!! In addition, the onset temperature for 
strongly anisotropic transport is comparable to the criti- 
cal temperature expected£3 and observedE-3 in other Ising 
quantum Hall ferromagnets. For future work it would be 
interesting to investigate transport in systems with de- 
liberately induced surface features. 
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